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Abstract 

We present a comprehensive structural characterization of two different highly 
pure nuclear graphites that compasses all relevant length scales from nanome¬ 
ters to sub-mm. This has been achieved by combining several experiments 
and neutron techniques: Small Angle Neutron Scattering (SANS), high- 
resolution Spin Echo SANS (SESANS) and neutron imaging. In this way 
it is possible to probe an extraordinary broad range of 6 orders of magnitude 
in length from microscopic to macroscopic length scales. The results reveal a 
fractal structure that extends from ~ 0.6 nm to 0.6 mm and has surface and 
mass fractal dimensions both very close to 2.5, a value found for percolating 
clusters and fractured ranked surfaces in 3D. 

1. Introduction 

Graphite has been used as a neutron moderator in several types of nu¬ 
clear reactors from the Ghicago Pile 1 in 1942 to the more recent Very High 
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Temperature Reactor (VHTR) and High Temperature Gas-cooled Reactors 
(HTGR). This synthetic polygranular material has a very high chemical pu¬ 
rity and a complex microstructure, which affects the mechanical properties 
under extreme conditions and irradiation damage [T]. 

5 The crystallite structure and disorder of graphite at the atomic level can be 
investigated by neutron or X-ray diffraction [2], and the microstructure by 
TEM, SEM or optical microscopy [3H7]. On the other hand, the bulk meso¬ 
scopic structure of the pores can be explored by small angle scattering of 
X-rays (SAXS) or neutrons (SANS) [8]. Very first SANS measurements on 
10 non-irradiated and irradiated nuclear graphites were performed in the 1960s 
[9] and 1970s [TU1 - IT2] . These results have been reinterpreted recently [13] 
to disclose a surface fractal structure from ~ 0.2 to 300 nm, i.e. over three 
orders of magnitude in length. However, the graphite inhomogeneities can be 
seen with an optical microscope or even with naked eye. Therefore an explo- 
15 ration over a larger range of length scales is necessary and for this purpose 
we have combined three neutron-based techniques: SANS, Spin Echo SANS 
(SESANS) and imaging to cover lengths from nm to mm. We investigated 
two different highly pure nuclear graphites, and the results show a fractal 
structure over an extraordinary large scale of lengths that spans 6 orders of 
20 magnitude and has fractal dimensions close to 2.5. This value is expected 
for several cases of percolating clusters [T1[I3] and in the most general case 
of fractured ranked surfaces [131 in three dimensions. 
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2. Methods 


2.1. Sample 

The samples were disk-shape specimens with a thickness of 0.5 mm and 
a diameter of 16 mm cnt from two types of nnclear graphite, designated 
5 as RID and PGA. The RID graphite was mannfactnred by Pechiney SA 
in the 1960’s by baking a paste made of oil coke and pitch, graphitized 
by electrical heating, and was nsed at the research reactor of the Reactor 
Institnte Delft. The PGA (Pile Grade A) graphite was mannfactnred by 
British Acheson Electrodes, Ltd. and Anglo Great Lakes, from needle shaped 
10 coke particles derived from the petrolenm indnstry. It was nsed in the early 
gas-cooled reactors in UK and has been object of several investigations [U 
m E]- PGA was mannfactnred by extrnsion, which leads to aligned coke 
particles along a direction e and thns to the anisotropic properties. For this 
reason we prodnced two series of samples: PGAl cnt perpendicnlar to e, 
15 and thus isotropic; and PGA2 cut along e and thus anisotropic. For the 
sake of simplicity in the following we will focus on PGAl and the results 
from PGA2 are averaged over the whole sample as for the other two samples. 
In this work, the supplement gives a detailed analysis of the results from 
the PGA2 sample and the anisotropy effect, which is weak but nevertheless 
20 visible. 

2.2. Neutron-matter interaction 

Since neutrons are electrically neutral, they can penetrate into matter 
deeply, only interacting with the nucleii, and are valuable probes of the struc¬ 
ture in the bulk. When a beam of thermal neutrons interacts with a sample 
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it may be scattered or absorbed with a respective probability that is given 
by the scattering lengths and the absorption cross sections of the specihc el¬ 
ements and isotopes in the sample (for a comprehensive introduction to neu¬ 
tron scattering see [SO])- These transmitted and scattered neutrons deliver 
5 structural information from the sample. In the following we will introduce 
the techniques used in our work. 

2.3. Small Angle Neutron Scattering (SANS) 

Fig. shows a schematic diagram of a SANS experiment. A neutron 
beam with a wavelength A is scattered by the structural inhomogeneities 
(pores/carbon matrix in the case of graphite) and recorded by a 2D position- 
sensitive neutron detector. The transmitted beam is captured by a beam-stop 
made of neutron-absorbing material, illustrated by the black circular area in 
the resulting scattering pattern in Fig. 

SANS measures the intensity of scattered beam I{Q), i.e. the scattering 
cross section, as a function of the scattering vector ^ that is related to the 
scattering angle 6 through Q = \ ^\ = 47rsin(6'/2)/A. I{Q) can be factorized 

as PU] : 

I{Q) = B-P{Q)S{Q). (1) 

Here H is a pre-factor given by the neutron scattering length density contrast, 
in our case between the carbon matrix and the pores. P{Q) is the form factor 
10 characterizing the morphology/shape of the pores, and S{Q) is the structure 
factor corresponding to the correlations between pores. 

In this work SANS measurements were performed at two instruments, the 
medium resolution PAXE and the high resolution TPA |2I] of the Laboratoire 
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Figure 1: (a) Schematic diagram of a SANS experiment. A collimated monochromatic 
(wavelength A) neutron beam is scattered by structural inhomogeneities in the sample 
and recorded by a 2D position-sensitive neutron detector. The magnitude of the scattering 
vector is Q = 1^1 = 47r sin(0/2)/A, with A the neutron wavelength, and 9 the scattering 
angle, (b) Schematic diagram of a SESANS experiment. A polarized monochromatic 
neutron beam passes through the setup, where two magnetic helds in opposite direction 
before and after the sample induce Larmor preccessions. The initial polarization state is 
completely recovered for the transmitted beam as illustrated by the blue path, whereas 
scattering changes the final polarization state of the beam, (c) Schematic diagram of a 
neutron imaging experiment. An incident polychromatic neutron beam is attenuated by 
the sample and the image recorded by a 2D camera-type neutron detector visualizes the 
structural inhomogeneities. 
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Leon Brillouin (LLB), CEA Saclay, France. On both instruments the neutron 
beam had a monochromatization of AA/A = 10%. The experiments on TPA 
were done at A=0.6 nm, covering the Q range 6 x 10“^ < Q < 1 x 10“^ nm“^, 
and on PAXE at A=0.37, 0.6, and 1.7 nm, respectively, covering the Q range 
5 3 X 10“^ < Q < 5 nm“^. The PAXE data for A=0.6 nm were brought to 
absolute units by normalization to the incident beam and were then used to 
normalise all other data using the large overlap in the Q-ranges illustrated 
by Fig. 1^ for the case of PGAl. 

2.4- Spin Echo Small Angle Neutron Scattering (SESANS) 

The basic principles of SESANS can be found in [22], and Fig. rep¬ 
resents schematically this technique. In contrast to conventional SANS, 
SESANS reaches high resolution accessing structural information on mi¬ 
crometres length scales by using a polarized neutron beam. Larmor pre¬ 
cessions are induced in two regions with opposite magnetic helds, and if the 
setup is symmetric, the beam recovers its initial state, leading to a maximum 
spin echo polarization. Scattering from the sample breaks the symmetry and 
reduces the echo polarization. SESANS measures the spin echo polarization 
Ps{z) as a function of the spin echo length z [23], the direction of which 
is determined by the geometry of the setup and it is always perpendicular 
to the propagation direction of the neutron beam n (see Fig. 03) • Ps{z) 
measures the projected scattering length density correlation function G'{z)\ 

Ps{z) = exp[fA^(G'(^) - G'(0))/27r], (2) 
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and G'(z) is the Hankel transformation of the SANS cross section: 

poo 

G'{z) = / MQz)I{Q)QdQ, (3) 

Jo 

where A is the wavelength of the nentron beam, t is the sample thickness, 
and Jo is a zeroth-order Bessel fnnction of the hrst kind. 

The high-resolntion nentron Spin Echo SANS experiments were performed 
on the dedicated instrnment of the Reactor Institnte Delft [22] at A=0.205 
5 nm, with AA/A = 5% and covered length scales from 30 nm to 20 /rm. 

2.5. Neutron imaging 

Fig. P shows a schematic view of a nentron imaging experiment. The 
incident polychromatic nentron beam passes throngh a sample. The inten¬ 
sity is attennated dne to absorption and/or scattering and is recorded by a 
10 2D camera-type nentron detector. Since different phases in the sample have 
different attennation coefficients, the resnlting image visnalizes the strnctnral 
inhomogeneities of the sample. 

Nentron imaging was performed on the cold nentron facility, ICON, of the 
Panl Scherrer Institnte, Switzerland PI- The samples were placed as close 
15 as possible to the detector and a resolntion of ~ 30 /xm was obtained. The 
transmission images were normalized to the empty beam and treated to ob¬ 
tain 8-bit digital images. These were farther analysed with the DIPimage 
Matlab toolbox (http://www.diphb.org/) as described in the snpplement. 
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Figure 2: (a) Absolute scattering cross section I{Q) of PGAl for the indicated different 
instrument and neutron wavelengths. The power law I{Q) oc Q~^ with /3 = 3.45 is 
illustrated by the black line, (b)-(d) Contrast variation experiments to investigate the 
effect of multiple scattering: (b) I{Q) of PGA2 in air (red circles) and of the same sample 
embedded in deuterated toluene (black triangles). The data are vertically shifted to overlap 
in (c). A slight difference is found at the very low-Q range, which is enlarged in (d). 










3. Results 


3.1. SANS 

The absolute scattering cross section of PGAl, shown in Fig. |^, fol¬ 
lows a power law I{Q) oc Q~^ with f3 = 3.45 ± 0.01. Similar results have 
5 been obtained for the other samples and the values of /3 are given in Table 
Deviations at high-Q’s are mainly due to (spin-incoherent) background, 
whereas at low Q’s to multiple scattering. Graphite is a strong neutron 
scatterer, a property used for nuclear applications, and multiple scattering 
is not negligible even for the thin samples used in this study, as illustrated 
10 by the transmission valnes given in the supplement. For this reason SANS 
was also measured by embedding one sample (PGA2) in denterated toluene, 
which reduces the scattering contrast. As a consequence the transmission 
increased and the scattered nentron intensity decreased by almost an order 
of magnitnde as shown in Fig. [^. However, the shape stays nnchanged with 
15 the exception of the very low-Q range illustrated by Fig. B. where the two 
curves are shifted on the log-log scale. In the enlarged view of Fig. |2]i, it is 
clear that the bending characteristic for mnltiple scattering [25] is less pro¬ 
nounced for the sample in deuterated toluene. Thus multiple scattering has 
an influence on the SANS patterns but does not affect the power law, which 
20 is the same for both cases and reflects the genuine strnctnral properties of 
the sample. Similar behaviour has been fonnd in expanded graphite [T8| and 
sedimentary rocks samples [2S]- 


9 




Figure 3: Normalized SESANS polarization Ps{z) of (a) PGAl and (b) RID graphite 
samples for different sample thicknesses (0.5, 1.0 and 2.4 mm) as a function of the spin 
echo length z. This effect, which due to the scattering power of the sample and thus to 
multiple scattering, is taken into account in (c), where Ps{tY°^* is plotted versus z, for 
Iq = 1 mm revealing the same generic curve for all samples. The solid lines in (a) and (b) 
correspond to the fitting curves from Hankel transform of Eq. and with the parameters 
of Table [21 
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3.2. SESANS 

Fig. 1^ shows the SESANS pattern of PGAl. As already mentioned, 
SESANS measures the Hankel transformation of I{Q) in the form of nor¬ 
malised (to the direct beam) SESANS polarization Ps, which is a function 
5 of the spin echo length z [23] • It probes length ranges from ~ 30 nm to ~ 20 
/im, and thus ’’sees” very large objects, which scatter a large fraction of the 
incoming beam. This is also the case in Fig. |^, where even for the 0.5 mm 
PGAl sample at large spin echo lengths z, the entire beam is scattered and 
Psi z > 10 nm) —)■ 0. However, even in this extreme case multiple scattering 
10 does not alter the results. This effect can indeed be accounted in a way simi¬ 
lar to the Beer-Lambert’s law in optics and if Ps{to) is the SESANS signal for 
a reference sample thickness to, the SESANS signal for any other thickness 
t is given by: Ps{ty°^* = Psito) [27|. This property was tested on three RID 
samples with thicknesses of 0.5, 1 and 2.4 mm respectively. As shown in 
15 Fig. 1^, Ps{z) decreases much faster for the thicker samples. However, when 
plotting Ps{ty°^^, with tQ=l mm, all data follow on the same generic curve 
shown in Fig. |^, with deviations when Ps is smaller than the experimental 
error of 10“^. More importantly, these results show that the porous structure 
extends up to macroscopic length scales, which can be investigated by direct 
20 imaging. 

3.3. Neutron imaging 

Fig. 1^ shows the neutron transmission image of the PGAl sample (the 
images for the other samples are in the supplement). The treatment of 
this image (see supplement) leads to a binary image, where the pores are 
25 distinguishable from the carbon matrix. Gonsequently the pores within the 
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Figure 4: : Neutron transmission image of the PGAl sample (a), within the selected yellow 
square, the pores are labelled with different colours (b), and a selected region is enlarged 
in (c). The analysis of the shapes of the pores leads to the Perimeter-Area plot of (d). 
The red line corresponds to a power law P cx with 7 = 0.77. The pore size distribution 
is given in (e), which shows a power law M{A) oc A'^ with r = —1.82 illustrated by the 
red line. 
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selected yellow square (of size 800 x 800 pixel) are labelled with different 
colours according to their sizes and shapes. This leads to Fig. i). and a 
selected region enlarged in (c) reveals a strongly ramihed, fractal structure. 
This structure may be brought in relation with the reported micropores 
s urn, that result from calcination and gas evolution during the manufacturing 
process and have sizes from several micrometres to hundreds of micrometres 
without preferred orientations. The binary images also lead to an estimation 
of the porosity 0image , which is about 82-85% of the porosity 0^) determined 
from density measurements, considering 2.25 g/cm^ as the density of single 
10 crystalline graphite. Thus pores smaller than the 30 pm (resolution limit of 
the images) contribute only by 15-18% to the total porosity. 

The plot of the pore perimeter P against area A for PGAl (the plots 
for other samples are given in the supplement), shown in Fig. |^, reveals a 
power law P (x. A"< with 7 = 0.77 ± 0.01 (red line in the hgure). This value 
15 is higher than 0.5, the value expected for Euclidian geometry, and is thus a 
signature of fractality. We note that in this hgure the points have different 
multiplicities. The pore size multiplicity M{A) is plotted against A in Fig. 
1^, where M{A) oc A'^ with r = —1.82 ± 0.03, illustrated by the red line 
in the hgure, and the values of r for other samples are listed in Table 
20 Thus the total volume of pores V{A) with specihc area size A is given by: 
V{A) oc M{A) X %3/2 oc A- X %3/2 = 7 lr+3/2 _ ^-0.3_ jpg gf ^pg 

exponent indicates that the pores with diherent sizes have almost the same 
total volume. 

A quantitative analysis of the transmission images involves a 2D-Fourier 
25 transform of the region inside the yellow square of Fig. which leads to 
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(a) PGAl 



Figure 5: :The 2D scattering pattern of PGAl obtained by 2D-Fourier transform on the 
selected region illustrated by the yellow square in Figj^. This leads to the I{Q) curve of 
(b). The power law at Q > 1.5 x 10“"^ nm“^, I{Q) oc illustrated by the green line, 

has an exponent smaller than 3 indicating a mass fractal. The red line in (b) corresponds 
to Eq. and with parameters of Table 

the 2D scattering pattern of Fig. 1^. By radially averaging this intensity, 
following the same procedure as for the SANS data, the scattering curve of 
Fig. 1^ is obtained, with Q calculated as given in the supplement. The 
resulting I{Q) levels off at low-Q region and crosses over to a power law 
5 I{Q) oc Q-2-44±0.09 pigPgj. 

4. Discussion 

These results disclose power laws, which are characteristics of fractal 
topology: fractal surfaces (of pores or particles) with a surface fractal di¬ 
mension Dg or mass fractals with a mass fractal dimension Dm [281129] . The 
most direct signature of fractality is the power law of the SANS intensity 
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Table 1: Power law exponents (/3, r), fractal dimensions {Dg, D^) and porosity and 
(/(image) derived from SANS and imaging. 



SANS 


Imaging 

Porosity 

Sample 

/? 

Dg 

Ds 

Dm ^ 

4 >d 

0image 

PGAl 

3.45 (1) 

2.55 (1) 

2.54 (2) 

2.44 (9) -1.82 (3) 

22.7(5)% 

18.6% 

PGA2 

3.43 (1) 

2.57 (1) 

2.50 (2) 

2.54 (9) -1.80 (4) 

22.7(5)% 

17.8% 

RID 

3.45 (1) 

2.55 (1) 

2.56 (2) 

2.45 (9) -1.74 (3) 

27.1(5)% 

23.0% 


shown in Fig. It is indeed expected that [ 301 1 ^ : 

for a snrface fractal: I{Q) oc (4) 

for a mass fractal: I{Q) cc (5) 

where the valnes of both and Dg are smaller than 3, the dimensionality 
of the Euclidian space. Expressing the power laws as I{Q) oc Q~^, (3 > 3 
corresponds to a surface fractal and /3 < 3 to a mass fractal. All graphites 
showed f3 > 3 leading to the surface fractal dimensions listed in Table 
Additional conhrmation of these results comes from the Perimeter-Area 
power law P oc of Fig. 111. where the exponent 7 is directly related to 
the surface fractal dimension Dg through 7 = {Dg — l)/2 [321I34] . leading to 
Dg = 2.54 ± 0.02 for PGAl. Similar results were obtained for the samples 
and all deduced values of Dg are listed in Table and are consistent with 
those obtained from SANS. Therefore imaging and SANS ’’see” the same 
surface self-similar fractal structure although they probe length scales more 
than 3 orders of magnitude apart. 

On the other hand, the power law of I{Q) derived from the Fourier trans- 
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formed images leads to an exponent significantly smaller than 3 (Fig. if), 
which is in line with a mass fractal and not snrprising at this very low Q- 
limit. The dednced valnes of Dm are given in Table [T| and for all samples 
D ^ D 

The experimental resnlts indicate that the pores of the graphite samples have 
a fractal (i.e. rongh) snrface while their assembly forms a mass fractal. To 
be specihc, one can consider that the pore strnctnre consists of many pore- 
bnilding blocks, which characterize the network of the pore clusters with a 
mass fractal property at the length scale above the size of the primary block. 
Moreover, the pore building block itself is bounded by a rough surface with 
fractal morphology. In this case the scattering function reflects both mass 
and surface fractal properties and Eq. can be factorized as [311 [35]: 

= <;/)£) Ap^ dvr £^/3 (6a) 

F(Q) = (1 + Q2^2)(^»-6)/2 (6b) 

■ sin[(T)m - 1) arctan(Q.^)], 

where Ap (=7.5 ■ 10^°cm“^ for carbon-air) is the scattering length density 
contrast, Vp the volume of a primary pore building block, £ the associated 
length, P{Q) the form factor, characteristic of the surface fractal morphol¬ 
ogy of the pores, S{Q) the structure factor corresponding to the mass fractal 
5 structure and T the gamma function. Besides the fractal dimensions Dm and 
Ds, two characteristic lengths are introduced: an upper cut-off length ^ and 
the length of the primary pore-building block £, which marks the cross-over 
between mass and surface fractal scattering. 
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Table 2: Parameters for the fractal model assuming Dg = Dm- 


SANS Imaging SESANS Global 


Sample 

Ds = 

(fixed) 

S (^m) 

(fitted) 

1 1 

^ A 

B (cm b 
(fitted) 

S{Q = 0) 

(calculated) 

B (cm b 
(fitted) 

B (cm b 
(caleulated) 

PGAl 

2.55 

18.9 (4) 

2.27 (10) 

6.6(1) X 10“ 

134(15) 

7.56(5) X 10“ 

6.3(9) X 10“ 

PGA2 

2.57 

17.0 (5) 

2.64 (25) 

8.4(1) X 10“ 

74(5) 

1.87(1) X 10“ 

9.8(5) X 10“ 

RID 

2.55 

18.7 (4) 

1.87 (6) 

3.0(1) X 10“ 

214(17) 

4.64(4) X 10“ 

4.2(5) X 10“ 


In order to fit all experimental results with this model in the most reliable 
way we adopted the following strategy : (1) the values of Dg were fixed to 
those of Tablederived from SANS; (2) the values of D^n and ^ were derived 
from the scattering patterns of the Fourier transformed images. These are 
5 given in Table [T| which shows that Dm ~ Dg. For the sake of simplicity we 
assumed Dm = Dg in the following. (3) With Dm, Dg and ^ hxed, i was 
determined by htting the SESANS data, that £11 the gap between SANS and 
imaging and probe the cross-over between mass and surface fractal scatter¬ 
ing. For this purpose the Hankel transform was performed numerically on 
10 Eq. (combined with Eq. through Eq. then the data were fitted using 
Eq. [^with only two floating parameters: length i and the prefactor B, and 
the resulting values are given in Table The fits are illustrated by the solid 
lines in Fig. [^and describe excellently the experimental findings. 

The final step is to combine SANS, SESANS and imaging, which is done 
15 in Fig. 1^ In this figure the lines correspond to the best fit of Eq. [^with 
the parameters of Table The Q values for which = 1 and Qi = 1 
are also indicated. We stress that the model provides excellent quantitative 
description of the experimental findings. Indeed the prefactors B given in 
Table which are deduced from fitting all data (global) are in excellent 
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Figure 6: Scattering curves for (a) PGAl (b) PGA2 and (c) RID over for the whole Q 
range investigated. The circles represent the SANS data; the orange crosses are from the 
Fourier transformation of the neutron images; and the gray shaded zone represent the 
Q-range covered by SESANS. The black lines correspond to the best fit of the Eq. [^with 
the parameters of Table The dotted lines correspond to Q = 1/^ and Xji respectively. 


agreement with the calculated values from the fitted values of I through Eq. 


6 a Therefore, a consistent picture of the fractal microstructure of graphite 


is obtained over an extraordinary large range of length scales of 6 orders of 
magnitude (~0.6 nm < 2t:/Q < 0.6 mm). Table shows that the fractal 
5 dimensions and cut-off lengths ^ are almost the same for all samples, which 
might be because both PGA and RID are made out of petroleum coke dU. A 
trend is only seen in the values of £, which is the parameter that reflects the 
particularities of the microstructure and also the anisotropy of PGA2 (see 
supplement). It has been shown that the fractal anisotropy can be found in 
10 graphite at low densities (expanded graphite) and that in fact under densi- 
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fication such anisotropy is lost [THl US] . 

In the literature besides the micropores mentioned above, so-called Mro- 
zowski cracks have been reported lasiE], that result from the anisotropic 
thermal shrinkage of the layered graphite structure during cooling from graphi- 
5 tization temperatures. These have lengths from nm to pm that are in the 
range probed by SESANS and SANS. It is therefore tempting to attribute 
the surface fractal scattering to these cracks, as suggested recently [13], and 
i to the cross-over between cracks and micropores. 

In contrast to optical or electron microscopy, the scattering techniques used 
10 in this work cannot distinguish between cracks and pores. However, they 
provide a statistical average of the correlations over the (macroscopic) sam¬ 
ples and reveal the most generic features of the structure. In this way it is 
possible to describe the complex and poly-disperse patterns of Fig. |^with the 
scattering law of Eq. involving a restricted number of parameters, which 
15 enables quantitative comparisons between different samples and systems. 
Besides the extraordinary large length scales over which fractality has been 
observed, the fractal dimensions found in this work are comparable to those 
expected for percolating clusters [HEB], for which heuristic arguments sug¬ 
gest that mass fractals should be bounded by their own natural fractal surface 
20 and Dm = Dg na. The porosity values of Table indicate indeed a topology 
close to that of a percolating cluster [3B] . 

Similar fractal exponents have been deduced for fracturing ranked surfaces 
in 3D [16], that could serve as a model for Mrozowski cracks. Therefore the 
identity Dm = Dg 2.5 is not coincidence but the consequence of the high 
25 degree of disorder, ramihcation and connectivity of the pore structure. Under 
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neutron irradiation surface fractality disappears |T3], and similar behaviour 
may be expected for the oxidised samples. The methodology developed in 
this work thus can be applied to further investigate the effect of irradiation 
damage and/or oxidation on the structural properties of graphite. 

5 Fractal scattering has also been reported for carbon nanopores [371 EH] , rocks 
[26] or cement [39] . The particularity of this work is in the extraordinarily 
broad length scale of six orders of magnitude over which fractal scaling is 
quantitatively valid. The combination of several techniques, from imaging to 
scattering and the methods can be applied to the investigation of other com- 
10 plex systems with a hierarchy of length scales such as biological materials, 
concrete and rocks, materials for C02 sequestration, Li batteries, fuel cells 
or solar cells. 
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